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ABSTRACT 

An algorithm of a convectional factory tray dryer with performance evaluation on the functionality of 
the dryer was reported. This is required to reduce drying time of material as well as man-machine interface 
which can cause slow operation and contamination of the material been dried. The algorithm of the model was 
developed with an interface accorded with computer software (MATLAB) for its quick implementation. This 
algorithm was able to predict parameter values needed for dryer automation such as heating element diameter, 
heating chamber length, pumping power and flow rate of air. Result generated from the algorithm, was found to 
be of importance in fabrication of a convectional factory electric tray dryer. The algorithm was also validated 
by manually calculated parameter values. The results were found to be the same, this proved the effectiveness of 
the algorithm. 
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I. INTRODUCTION 

The term "drying" is basically utilised in thermal engineering environment which may include food or 
wood [1]. Drying is an ancient means of preserving food [2]. It is so important in industry because the kinetics 
of the process can define the final quality properties of the dried material. Haghi and Amanifard, [3] reported 
that it is a complex process involving simultaneous coupled transient heat, mass and momentum transport. In 
domestic and industrial parlance, there are many types of dryers, such as: drum dryers, tunnel dryers, tray 
dryers, drum dryers, fluidized bed dryers, spray dryers, flash dryers, belt dryers, rotary dryers, freeze dryers, and 
vacuum dryers [2]. Of these dryers, tray dryer is more economical and not complex. The materials to be dried 
are easily and conveniently dispersed along the surface area of the tray. Drying operations can be broadly 
classified according to whether they are batch or continuous [4]. These terms are applied specifically from the 
point of view of the substance being dried. The process term batch drying entails the quantity of the substance 
being dried passing through a continuous stream of hot flowing air. In the continuous process the substance to 
be dried as well as the gas passes continually through the equipment. The mechanism involved in batch drying 
include; constant rate period, falling rate period and equilibrium moisture content. According to Shepherd [5], 
the rate of evaporation and the surface temperature can be obtained by heat balance. Foust et al.,[6] reported on 
the graphical models that can be used to describe drying operations. Akpinar, [7] concluded that at every mode 
of drying route under exanimation, there was no constant-rate period of drying. The rate of drying was reported 
to be decrease along with the drying time [8]. Economical and rational uses of agricultural produce demands a 
means for the conservation of excess production. Thus, the need for an effective drier to dry to a predetermined 
level the moisture contained in agricultural produce such as meat, fish cassava, potatoes. Hence, the present 
work seeks to establish an algorithm for a convectional factory tray dryer. 

1.1 DESIGN CONSIDERATIONS AND ANALYSIS 

The automated tray dryer was designed to have the following nomenclature: heating element, the prime 
mover, fan (blower), control system and algorithm (flow chart). 
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1.1.1 Overview of Heating Chamber 



Pi,Tipi,Ui 



Control Volume 
(heating chamber) 



T 



P2, T2, p2, U2 



q =Cp (To2- Toi) 

Design of Heating Element 

Heat is radiated from the coil and by means of forced convection is transported to the commodity 
through a blower. Coil of circular-cross sectional or rectangular ribbons was use as heating elements. Under 
steady state condition, a heating element dissipates nearly as much heat from its surface as it receives power 
from the electric supply Theraja,[9]. If P if the power input and it is the heat dissipated by radiation, hence, 
P =H (1) (under steady state condition) 

From Stefan's law of radiation, heat radiated by a hot body is given by; 
H=5 .72e^ [(Ti/100)^- (T2/100)'^] wW 

Now P = ^ and R = p(L/A) = ^ . Therefore, P = Tid^ vV (4pZ) or Z/d^ = nx^l{ApP) 
Where v = the potential gradient (in volts), R= the resistance of the heating element. 

Total surface area of the coil element = udl (2) 
If H is the heat dissipated by radiation per seconds per unit area of the coil, the heat radiated per second is; 

P = (ndlH ) or ^^^-^ = nd^lH 

(4pO 

Where: d = diameter, l=length, p = density and H = heat dissipated 

Therefore; — = — — 
I 

(3) 

1.1.2 Fan Motor Relationship 

The fan modeled for this purpose is an axial fan. The model is a function of the following parameters; 
Density of fan air 

Fan air density is a mass per unit volume of the air. The density of a perfect gas is a function of its molecular 
weight, temperature and pressure as indicated by 

= (M/386.7)(529.7/r,)(&;,/29.92) (4) 

Where M = molecular weight, T^= temperature and bx = pressure 
For dry air, the above equation reduces to 

p^ = 132S{bjT^) (5) 

This equation is valid for moist air. 

Capacity of fan 

Since the fan is of the axial type, most of our design parameters would be limited to analysis in x- direction 

therefore; = 1.097i4^V(P^x/Px) (6) 
Thus the total capacity, which is the volumetric flow rate at fan air density Py-, 

is given as Q = QxPx/P (7) 

A more practical value of Q is given as Q — CodR (8) 
The value of Q is limited as demonstrated; A blade has a mean sectional area A and length CcR thus, the volume 
of blade becomes OCRA, its acceleration; (x) R and 

Centrifugal force F =p{aRA)o)^R stress due to force is F / A =0C pco^R^ (9) 

If Sf is the failure stress blade material then the content 5 < Sj^/u (10) 

Combining equations (9) and (10), we have 

1 1 1 

0) < [Sf/anpY V/? and Q < C/? V(<^n)2 [5^/p]^ (11) 
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Fan total pressure 

Fan total pressure is the differences between the total pressure at the fan outlet and the total pressure at the fan 
inlet. 

Pt = Pt2 - Ptl (12) 
If the fan draws air directly from the atmosphere as in this case 

Ptl = 0, i.e Pt = Pt2 (13) 
If also, the fan is connected to a duct as in the case of our model, the pressure drop AP2-X 
must be added i.e + AP2-X (14) 

Fan velocity pressure/ Fan static pressure 

Fan velocity pressure is the pressure corresponding to the average velocity at fan outlet 

Pv = (Q2/109.7/l2)2p2 (15) 

The Fan static pressure; this is the difference between the fan total pressure and the fan velocity pressure. 
Therefore, fan static pressure is the difference between the static pressure at the fan outlet and total pressure at 

fan outlet. = - = P^^_P^^ (16) 
If the fan draws air directly from the atmosphere P^ = P^ — Py = P^^ 

Fan speed 

The fan speed simply refers to the relative speed of the impeller. 
Compressibility factor 

This is the ratio of the fan total pressure P^ that would develop for an incompressible fluid to the fan total 
pressure P^ that is developed for a compressible fluid when all other conditions remain the same. 

= P't/Pt = n(n - l)[(Pt2a/^tla)""'^" - '^]/iPt2a/Ptla) (17) 
Compressibility factor may also be determined by measurement using 

X = Px/Ptla (18) 

Z = [(/ - l)/y] 6356W(QPna) (19) 

Then = [zlogil + x)]/[x log (1 + z)] (20) 

Fan power input, output / total fan efficiency (jJt) 

This is the power require to drive the fan and elements in the drive train which is considered a part of the fan. 
Fan power input can be calculated from appropriate measurements using a dynamometer, torque meter or 
calibrated motor. 

However, the fan power output is expressed as Hq = QP^K^/63S6 (21) 
The total fan efficiency is given as; 

3r],= QP,K^/(63S6H) (22) 
or rjt = a)J a)sn 

Fan static efficiency tj^ 

The fan static efficiency is expressed as a function of the total fan efficiency, fan static pressure and fan total 
pressure. The relationship is expressed below; 

Vs = VtPs/Pt (23) 
Fan sound Level 

The sound Level is ten times the logarithm (in base ten) of the ratio of the actual sound power in watt to 
10"^2 watts. Ly, = 10\og(w/10-^^) (24) 



1.1.3 Carrying capacity of tray 

The slot carrying the frame will be welded to the main frame of the drying chamber. The characteristic of the 
welded tray supports and of the tray itself are given below. 
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Unit of polar 


A= 1.414bd 


X = b/2 


moment of inertia 




Y = d/2 


\ =d(3bVd^)/6 



Figure 1 . (weld diagram and design considerations of the tray section) 

Primary shear force produced in the weld t' = Sf/ A (25) 
The Moment Mi due to the weight of the commodity to be dried (assume to be concentrated at the centre of the 

plate) produces a normal bending stress Tjy in the welds. It is normal to assume that this stress acts normal to the 
throat area, during weld analysis. 

Unit moment of inertia = d(3b^ + d^)/6 (26) 
Moment of inertia on throat area 

Jt = 0.707/^ (27) 

S=Mc/l (28) 

The factor of safety of guarding against static failure in the weld metal is Tl=Sy/ S (29) 

1.1.4 The Control System 

The control system is an integral system that monitors the condition within the drying section of the 
machine such that when the condition indicates a dry state, it by means of manipulation of a programmable logic 
and control (PLC) network, sets off the shutdown process. The control system is a logic assembly of the major 
components of the machine such as the actuator motor and the heating element. The sensor which is an 
electronic device monitors the moisture level within the drying chamber and interprets it is as a voltage drop and 
then sends this signal to a comparator which is compared with the voltage drop across the motor coil when a 
zero voltage is recorded at the comparator over a time stretch, the time on De -energizing control units sends a 
confirmatory shutdown signal to the power unit. Meanwhile, a unit feedback loop continually relates the system 
response signal back to a summation point (or a comparator) which also works like the former. The control 
system, in effect, is a two way system which maximizes its performance. A moisture detector could be fitted 
into the drying chamber in order to measure the moisture level in any commodity that is been dried. When the 
commodity is sufficiently dried, the moisture detector will trigger a system (by sending an analogue or digital 
signal is a robotic arm) for instance, that will turn off the drier and the blower. The system arrangement is 
diagrammatically shown below: 



Linkv 



Moisture detector P»vot 
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K2 
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Figure 2. Control system arrangement 
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1.1.5 Heating Element Dimensioning 

Nichrome wire and Nickel-Chrome are most widely used heating element. They can be safely heated above 

1000°C and their emissivity is 0.09 each. Their radiating efficiency is in the neighbourhood of 0.5. For primary 

design calculations, the following holds: 

Heating element material: Nichrome 

Charge of the heating element: 30°C 

Final temperature: 800°C 

Radiating efficiency: 57% 

Q 

Specific resistance of Nichrome: 109x10" Qm 

Source of electricity supply (domestic use): 220V, 650W 

l/d^ = nv^/4pP (30) 
L/d = 53653170.23 

H = 5.72 ek[(TjlOOy - (72/100)^] = 38870 
H = 11.49m and d = 0.0045m = 4.5mm 



1.1.6 Heating Chamber Length and Pumping Power of Blower 

The cross- section of the heater tube is given as D = 200 mm = 0.2m 
The air temperature reaching the commodity was about 40°C 

The air velocity is about 3 — SlTlls 
The air entering the blower is at 15°C 
The Reynolds analogy will be followed in the design process 

The friction factor is calculated by f = 0.0791 
The air properties are considered at the mean temperature 

Atn = (Ati - At2)/ln(Ati/At2) 

II. RESULTS AND DISCUSSION 

Calculation: 

The mean film temperature is given as tj^ = (t^ + = 413.15°C 

The properties of air at tj^ = 686.75 (from table) 
We have prandtl number = 0.681 
Kinematics viscosity of air, fc j = 3.2 X 10"^ mVs 
Air density, p^j^ = 0.73kg /m^ 

Thus, Reynolds number Re — CD / Ki 

(33) 

Re = 31250 

i?g 31250 > 2000, Confirms that the flow is turbulent 

F = 0.00595 

Staton number is given as Sf = 

(34) 

s,= ^ = 0.002975 

2 

Nu = 63.312 

a = Nu*A/d = 0.04 

Mass flow rate of air m = pAV = 0.57kg /s 
Where Q = oc AAt 

(35) 

Q = mc{ta2 - tfli) 

(36) 

Q= 14.72/cw 
At„ = 772A: 
L = 72cm 
A = 1.70cm^ 



(31) 
(32) 
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The Pumping Power of the Fan / Q = / cAt 

^ = Ppumping = 11-60 W 
Flow rate of pumped air 

Ho = 

The SI unit system this equation reduces to P = QPf^Kp and 

Q = 0.34 mV^ 

Free body diagram and shear force diagram 

Figure 3 and 4 are representative of the free body and shear force diagram respectively. 



(37) 





Figure. 3. The free body diagram at tray supported 



h 



h 



The figure show the free body diagram at tray supported in both end by a supporting slot welded to the wall of 
the drying chamber. The free body diagram of one end of the tray will show a shear force Sf and reaction 
moment M (see below). 
V 



M 



The primary shear in the welds is of magnitude 
A= 1.414hb 

Chosen b = 45 cm = 450 mm and h = 5mm 
A = 3181.5 mm^ 

Thus, with a maximum weight of 10 kg equally supported by both end Sf = 49.05N 

The primary shear stress is t' = O.OISM PCL 
The bending moment diagram is shown below; 



(38) 
(39) 
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(40) 



The maximum bending occurring at the center of the tray is given as 
M= 9800 Nmm, 

The unit polar moment of inertia is 
^ _ d{3d^+d^) 

5u = 3042000 mm^ 

The polar moment of inertia base on the throat area (equation) is given as 

/ = 707hS^ (41) 
Substituting values of h and, / = 10753470 mm^ 

The centroid of the weld is located at 

X = b/2 = 225mm 
Y = d/2 = 15mm 

We can fine the secondary shear stress in components parallel is x and y, the component of y is; 

Ty = — (42) 

Ty = 0.2Mpa 

The component x is; 

ri' = MY /J = OmSGMPa 

These stress components should be combined to yield the maximum stresses; which occur at both end of the tray 
in the welds 

^ = V^x + (43) 

T = 0.21 MPa 

Chosen AWS electrode E60xx, the tensile strength^ 62 kpsi 
Yield strength =50 kpsi 

Elongation percentage =17-26% (Data based on the America Welding Society specification standard), if we 
base the failure of weld on yielding, we have, 

5y = 50 kpsi = 3445 MPa. The factor of safety guarding against static failure of the weld is provided 

by n = Sy/S = 1640.4 

CONCLUSION 

The value indicating total safety is somewhat large. This indicates that weld can support a load much 
greater than the 10 kg considered in the design calculations. However, the value of the maximum load that can 
be carried should be decided taking into account that the tray will bend after a certain limit of loading. This limit 
can be gotten theoretically but limiting constraint as the strength properties of the material available in the local 
market will decide on the maximum loads that can be carried. This system is a trade-off between efficiency and 
reliability. The result could be easily implemented, if any output parameter found is unacceptable by the 
program user, a new set of parameters can be computed into the program until an optimal design is obtained. 
The system occupies a total space of about 600 x 700 x 500 mm and weight about 15 kg. A total of 30 kg can 
be dried at a time. 
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Heating 
Element 



Appendix 

FLOW CHART 




State the values of; 
Xj Ij dj ti and 



Compute; 

Tp =tfi+ 273 



State values of; 

Ki= PGnematic viscosity of air at me an film temp. 
Pr = Prandtl number at mean film temp. 
C = Air velocity through the system 
D = Direction of heating chamber 



Determine Reynolds number (Re) 
Re = CB/Ki 



Heating 
Chamber 




Turbulentflow 
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Compute F, for turbulent floT;^^; 



Heating 
Chamber 



State values of; 

Thermal conductivity of fluid. 
p2 = Ivlass density at mean film temp. 
= Total lensth of heating chamber 





C o mpute ; 
d = A/D 

_ _ CTi - ti!) - CT3 - ta!) 




A = H^/dAt 

J^mp = C^/1.027At * 



State values of, 

Pt= Differential pressure across the bloT/^^er. 
Kp = Compressibility factor of the fan 



Flow 



C o mpute ; 
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